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Abstract 

A parking function of length n is a sequence {bi,b2, ■ ■ ■ ,bn) of nonnegative integers for which 
there is a permutation tt S S'„ so that < < i for aU i. A well-known result about 

parking functions is that the polynomial Pn{q), which enumerates the complements of parking 
functions by the sum of their terms, is the generating function for the number of connected 
graphs by the number of excess edges when evaluated at 1 + q. In this paper we extend this 
result to arbitrary connected graphs G. In general the polynomial that encodes information 
about subgraphs of G is the Tutte polynomial tQ{x, y), which is the generating function for two 
parameters, namely the internal and external activities, associated with the spanning trees of 
G. We define G- multiparking functions, which generalize the G-parking functions that Post- 
nikov and Shapiro introduced in the study of certain quotients of the polynomial ring. We 
construct a family of algorithmic bijections between the spanning forests of a graph G and the 
G-multiparking functions. In particular, the bijection induced by the breadth-first search leads 
to a new characterization of external activity, and hence a representation of Tutte polynomial 
by the reversed sum of G-multiparking functions. 
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^Partially support by NSF grant DMS 0245526. 



1 



1 Introduction 



The (classical) parking functions of length n are sequences (61, 62, . . . , 6„) of nonnegative integers for 
which there is a permutation vr € S'n so that < < i for all i. This notion was first introduced 
by Konheim and Weiss in the study of the linear probes of random hashing function. The 
name comes from a picturesque description in of the sequence of preferences of n drivers under 
certain parking rules. Parking functions have many interesting combinatorial properties. The most 
notable one is that the number of parking functions of length n is (n + Cayley's formula for 

the number of labeled trees on n + 1 vertices. This relation motivated much work in the early study 
of parking functions, in particular, combinatorial bijections between the set of parking functions of 
length n and labeled trees on n + 1 vertices. See |H1 for an extensive list of references. 

There are a number of generalizations of parking functions, for example, see jl] for the double 
parking functions, |2U1 122105] for /s-parking functions, and jl7lll4j for parking functions associated 
with an arbitrary vector. Recently, Postnikov and Shapiro ^Hj proposed a new generalization, the 
G-parking functions, associated to a general connected digraph D. Let G be a digraph on n + 1 
vertices indexed by integers from to n. A G-parking function is a function / from [n] to N, the 
set of non-negative integers, satisfying the following condition: for each subset U C [n] of vertices 
of G, there exists a vertex j ^ U such that the number of edges from j to vertices outside U 
is greater than /(j). For the complete graph G = Kn+i) such defined functions are exactly the 
classical parking functions, where one views i^n+i as the digraph with one directed edge (i, j) for 
each pair i 7^ j. In Pj Chebikin and Pylyavskyy constructed a family of bijections between the set 
of G-parking functions and the (oriented) spanning trees of that graph. 

Perhaps the most important statistic of the classical parking functions is the (reversed) sum, 
that is, (2) — (xi + X2 + • • • + Xn) for a parking function (xi, X2, . . . , Xn) of length n. It corresponds 
to the number of linear probes in hashing functions |lUj . the number of inversions in labeled trees 
on [n + 1] ^3], and the number of hyperplanes separating a given region from the base region in the 
extended Shi arrangements HOj , to list a few. It is also closely related to the number of connected 
graphs on [n + 1] with a fixed number of edges. In the second author gave a combinatorial 
explanation, which revealed the underlying correspondence between the classical parking functions 
and labeled, connected graphs. The main idea is to use breadth-first search to find a labeled tree 
on any given connected graph, and record such a search by a queue process. 

The objective of the present paper is to extend the result of to arbitrary graphs. For a 
general graph G, a suitable tool to study all subgraphs of G is the Tutte polynomial. This is a 
generating function with two parameters, the internal and external activities, which are functions 
on the collection of spanning trees of G. Evaluating the Tutte polynomial at various points provides 
combinatorial information about the graph, for example, the number of spanning trees, spanning 
forests, connected subgraphs, acyclic orientations, subdigraphs, etc. These many valuations make 
the Tutte polynomial one of the most fundamental tools in algebraic graph theory. 
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An important approach to getting information about the Tutte polynomial is to use partitions. 
This approach dates from the 1960's, see Crapo j^. More information on the history of Tutte 
polynomial can be found in 9 . Also in that paper Gessel and Sagan proposed a number of 
new notions of external activity, along with a new way to partition the substructures of a given 
graph. The basic method is to use depth-first search to associate a spanning forest F with each 
substructure to be counted. This process partitions the simplicial complex of all substructures 
(ordered by inclusion) into intervals, one for each F. Every interval turns out to be a Boolean 
algebra consisting of all ways to add external active edges to F. Expressing the Tutte polynomial 
in terms of sums over such intervals permits one to extract the necessary combinatorial information. 

In [U] Gessel and Sagan also mentioned another search, the neighbors-first search, and related the 
external activity determined by the neighbors-first search on a complete graph with n + 1 vertices to 
the sum of (classical) parking functions of length n. This connection was further explained in j23j . 
In the present paper we extend this result to an arbitrary graph G by developing the connection 
between Tutte polynomial of G and certain restricted functions defined on V{G), the vertex set of 
G. This is achieved by combining the two approaches mentioned before. First, we use breadth-first 
search to get a new partition of all spanning subgraphs of G. Each subgraph is associated with a 
spanning forest of G, which allows us to get a new expression of the Tutte polynomial in terms of 
breadth-first external activities of its spanning forests. Second, we construct bijections between the 
set of all spanning forests of G and the set of functions defined on V{G) with certain restrictions. 
One of such bijection, namely the one induced by breadth-first search with a queue, leads to the 
characterization of the (breadth-first) external activity of a spanning forest by the corresponding 
function. 

To work with spanning forests, we propose the notion of a G-multiparking function, a natural 
extension of the notion of a G-parking function. Let G denote a graph with a totally ordered vertex 
set V{G). Often we will take V{G) = [n] = {1,2,... ,77,}. For simplicity and clarity, we assume 
that G is a simple graph in most of the paper, except at the end of Section 3 where we explain 
how our construction could be modified to apply for general directed graphs, with possible loops 
and multiple edges. This includes undirected graphs as special cases, as an undirected graph can 
be viewed as a digraph where each edge {u,v} is replaced by a pair of arcs {u,v) and {v,u). 

For any subset U C V{G), and vertex v & U, define outdegu{v) to be the cardinality of the set 
{{v,w] G E{G)\w i U}. Here E{G) is the set of edges of G. 

Definition 1. Let G he a simple graph with V{G) = [n]. A G-multiparking function is a function 
f : V{G) = [n] ^ N U {00}, such that for every U C V{G) either (A) i is the vertex of smallest 
index in U, (written as i = ui\n{U)), and f{i) = 00, or (B) there exists a vertex i (z U such that 
< f{vi) < outdegu{i). 

The vertices which satisfy f{i) = 00 in (A) will be called roots of f and those that satisfy (B) 
(in U) are said to be well-behaved in U, and (A) and (B) will be used to refer, respectively, to 
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these conditions hereafter. Note that vertex 1 is always a root. The G-multiparking functions with 
only one root (which is necessarily vertex 1) are exactly the G-parking functions, as defined by 
Postnikov and Shapiro. 

Sections 2 and 3 are devoted to the combinatorial properties of G-multiparking functions. In §2 
we construct a family of algorithmic bijections between the set MVg of G-multiparking functions 
and the set Tg of spanning forests of G. Each bijection is a process based on a choice function, 
(c.f. §2), which determines how the algorithm proceeds. In §3 we give a number of examples to 
illustrate various forms of the bijection. This includes the cases where there is a special order on 
V{G), for instance depth-first search order, breadth-first search order, and a prefixed linear order, 
and the cases that the process possesses certain data structure, such as queue and stacks. At the 
end of §3 we explain how the algorithm works for general directed graphs. 

Section 4 is on the relation between G-multiparking functions and the Tutte polynomial of 
G. First, for each forest F we give a characterization of F-redundant edges, which are edges of 
G — F that are "irrelevant" in determining the corresponding G-multiparking function. Using that 
we classify the edges of G, and establish an equation between |-E(G)| and sum of G-multiparking 
function, \E{F)\, and the F-redundant edges. Then we use breadth-first search to partition all 
the subgraphs of G into intervals. Each interval consists of all graphs obtained by adding some 
breadth- first externally active edges to a spanning forest F. The set of breadth- first externally 
active edges of F are exactly the F-redundant edges of a certain type, which allows us to express 
the number of breadth-first externally active edges, and hence the Tutte polynomial, by the values 
of corresponding G-multiparking function. In section 5 we exhibit some enumerative results related 
to G-multiparking functions and substructures of graphs. 

2 Bijections between multiparking functions and spanning forests 

In this section, we construct bijections between the set A4Vg of G-multiparking functions and the 
set J-G of spanning forests of G. For simplicity, here we assume G is a simple graph with V{G) = [n] . 
A sub-forest F of G is a subgraph of G without cycles. A leaf of F is a vertex v € y{F) with degree 
1 in F. Denote the set of leaves of F by Leaf{F). Let Y\ be the set of all ordered pairs {F, W) 
such that F is a sub-forest of G, and 7^ C Leaf{F). A choice function 7 is a function from 
to V{G) such that 7(F, W) E W. Examples of various choice functions will be given in §3, where 
we also explain how the bijections work on a general directed graph, in which loops and multiple 
edges are allowed. As one can see, loopless undirected graphs can be viewed as special case there. 

Fix a choice function 7. Given a G-multiparking function / G MVg, we define an algorithm 
to find a spanning forest F € J^g- Explicitly, we define quadruples {vali, Pi,Qi, Fi) recursively for 
i = 0,1, . . . ,n, where vali : V{G) ^ Z is the value function. Pi is the set of processed vertices, Qi is 
the set of vertices to be processed, and Fi is a subforest of G with V{Fi) = Pj U Qi, Qi C Leaf{Fi) 
or Qi consists of an isolated vertex of Fi. 
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Algorithm A. 

• Step 1: initial condition. Let valQ = f , Pq he empty, and Fq = Qq = {1}. 

• Step 2: choose a new vertex v. At time i > 1, let f = 7(-Fi--i, where 7 is the 
choice function. 

• Step 3: process vertex v. For every vertex w adjacent to v and w ^ Pi-i, set vali{w) = 
vali-i{w) — 1. For any other vertex u, set vali{u) = vali-^i{u). Let N = {w\vali{w) = 
-1, vali^i{w) / -1}. Update Pi, Qi and Fi by letting P, = P,_i U {v}, Qi = Qi-i UN\{v} 
if Qi-i UN\{v} 7^ 0, otherwise Qi = {u} where u is the vertex of the lowest-index in [n] — Pi. 
Let Fi be a graph on Pi U Qi whose edges are obtained from those of by joining edges 
{w, v} for each w N. We say that the vertex v is processed at time i. 

Iterate steps 2-3 until i = n. We must have Pn = [n] and Q„, = 0. Define $ = $7,0 • ■M'Pg ~^ 
Tg by letting $(/) = F„. 

If an edge {u, vj\ is added to the forest Fi as described in Step 3, we say that w is found by 
V, and V is the parent of w, if v £ Pi-i- (In this paper, the parent of vertex v will be frequently 
denoted v^.) By Step 3, a vertex w is in Qi because either it is found by some v that has been 
processed, and {v, w} is the only edge of Fi that has w as an endpoint, or w is the lowest-index 
vertex in [n] — Pi and is an isolated vertex of Fi. Also, it is clear that each Fi is a forest, since 
every edge {u, w} in Fi \ Fi^i has one endpoint in V{Fi) \ V{Fi-i). Hence 7(-Fj, Qi) is well-defined 
and thus we have a well-defined map <I> from MVg to J-q- The following proposition describes the 
role played by the roots of a G-multiparking function /. 

Proposition 2.1. Let f he a G-multiparking function. Each tree component T of ^{f) has exactly 
one vertex v with f{v) = 00. In particular, v is the least vertex ofT. 

Proof. In the algorithm A the value for a root of / never changes, as 00 — 1 = 00. Each nonroot 
vertex it; of T is found by some other vertex v, and {v, w} is an edge of T. As any tree has one 
more vertex than its number of edges, it has exactly one vertex without a parent. By the definition 
of Algorithm A, this must be a root of /. 

To show that the root is the least vertex in each component, let ri < r2 < • • • < be the roots 
of / and suppose Ti,T2, . . . ,Tk are the trees of F = ^(/), where G Tj. Let Tj be the tree of 
smallest index j such that there is a i; € Tj with v < rj. Then j > 1 since the vertex 1 is always a 
root. Define U := ViTj UTj+i U . . . UT^). U is thus a proper subset of V{G) = [n]. By assumption, 
the vertex of least index in U is not a root. Therefore, U must contain a well-behaved vertex; that 
is, a vertex v such that < f{v) < outdegu{v). Note that all the edges counted by outdegu{v) lead 
to vertices in the trees Ti, T2, . . . , By the structure of algorithm A, all the vertices in the first 

j — f trees are processed before the parent of v is processed. But this means that by the time A 
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processes all the vertices in the first j — 1 trees, vali{v) = f{v) — outdegu{v) < —1, so v should be 
adjacent to some vertex in one of the first j — 1 trees. This is a contradiction. 

□ 

From the above proof we also see that the forest F = ^(/) is built tree by tree by the algorithm 
A. That is, if Tj and Tj are tree components of F with roots rj, Vj and < r^, then every vertex 
of tree Tj is processed before any vertex of Tj. 

To show that is a bijection, we define a new algorithm to find a G-multiparking function for 
any given spanning forest, and prove that it gives the inverse map of <1>. 

Let G be a graph on [n] with a spanning forest F. Let Ti,...,Tk be the trees of F with 
respective minimal vertices ri = 1 < r2 < • • • < r^. 
Algorithm B. 

• Step 1. Determine the process order vr. Define a permutation vr = (7r(l), 7r(2), . . . , vr(n)) = 
{viV2 ■ ■ ■ Vn) on the vertices of G as follows. First, vi = 1. Assuming vi,V2, ■ ■ ■ ,Vi are deter- 
mined, 

— Case (1) If there is no edge of F connecting vertices \nVi = {f i, • • • j Vi} to vertices 
outside Vi, let Vj+i be the vertex of smallest index not already in V^; 

— Case (2) Otherwise, let W = {v ^Vi : v is adjacent to some vertices in Vi}, and F' be 
the forest obtained by restricting FtoViU W . Let Uj+i = 7(F', W). 

(Hereafter, when discussing process orders, we will write Vi as vr(i).) 

• Step 2. Define a G-multiparking function / = fp. Set /(n) = f{r2) = ■ ■ ■ = f{rk) = 
oo. For any other vertex v, let be the minimal vertex in the tree containing v, and 
V, v'P, ui, . . . ,ut,ry be the unique path from v to r^. Set f{v) to be the cardinality of the set 

{v,\iv,v,)eE{G), 7T-\vj)<n-\vP)}. 

To verify that a function f = fp defined in this way is a G-multiparking function, we need the 
following lemma. 

Lemma 2.2. Let f : V{G) —>■ N U {oo} be a function. If v £ U CI V{G) obeys property (A) or 
property (B) and W is a subset of U containing v, then v obeys the same property in W. 

Proof. If f{v) = oo and v is the smallest vertex in U, then clearly it will still be the smallest 
vertex in W. If v is well-behaved in U, then < f{v) < outdegu{v) and as W Q U, we have 
outdegu{v) < outdegwiv)- Thus v is well-behaved in W. 

□ 
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The burning algorithm was developed by Dhar j7j to determine if a function on the vertex set 
of a graph had a property called recurrence. An equivalent description for G-parking functions is 
given in We mark vertices of G starting with the root 1. At each iteration of the algorithm, 
we mark all vertices v that have more marked neighbors than the value of the function at v. The 
function is a G-parking function if and only if all vertices are marked when this process terminates. 
Here we extend the burning algorithm to G-multiparking functions, and write it in a linear form. 

Proposition 2.3. A vertex function is a G -multiparking function if and only if there exists an 
ordering 7r(l), 7r(2), . . . , 7r(n) of the vertices of a graph G such that for every j, 7r(j) satisfies either 
condition (A) or condition (B) in Uj := {7r(j), . . . ,Tr{n)}. 

Proof. We say that the vertices can be "thrown out" in the order 7r(l), 7r(2), . . . , 7r(n) if they satisfy 
the condition described in the proposition. By the definition of G-multiparking function, it is clear 
that for a G-multiparking function, vertices can be thrown out in some order. 

Conversely, suppose that for a vertex function / : V{G) ^ N U {00} the vertices of G can 
be thrown out in a particular order '/r(l), vr(2), . . . , 7r(n). For any subset U of V{G), let k be the 
maximal index such that U Q Uk = {Tr{k),... ,7r(n)}. This implies TT{k) € U. But vr(A;) satisfies 
either condition (A) or condition (B) in Uk- By Lemma 12.21 Tr{k) satisfies either condition (A) or 
condition (B) in U. Since U is arbitrary, / is a G-multiparking function. 

□ 

Proposition 2.4. The Algorithm B, when applied to a spanning forest ofG, yields a G-multiparking 
function f = fp- 

Proof. Let vr be the permutation defined in Step 1 of Algorithm B. We show that the vertices can 
be thrown out in the order 7r(l), 7r(2), . . . , 7r(n). As 7r(l) = 1, the vertex 7r(l) clearly can be thrown 
out. Suppose vr(l), . . . ,7r(A; — 1) can be thrown out, and consider 7r(A;). 

If f{ir{k)) = 00, by Case (1) of step 1, 7r(A;) is the smallest vertex not in {7r(l), . . . ,7r{k — 1)}. 
Thus it can be thrown out. 

If f{'jr{k)) 7^ 00, there is an edge of the forest F connecting 7r{k) to a vertex w in {■7r(l), . . . , 7r{k — 
1)}. Suppose w = 7r{t) where t < k. By definition of /, there are exactly f{7r{k)) edges connecting 
7r{k) to the set {7r(l), . . . ,7r(t — 1)}. Hence f{7r{k)) < outdeg^T^(^k),...,TT(n)}{'^if^))- Thus TT{k) can be 
thrown out as well. 

By induction the vertices of G can be thrown out in the order vr(l), vr(2), . . . , 7r(n). 

□ 

Define ^^^^g : ^ MVg by letting ^'^,g(/) = fp- Now we show that — ^7,G' s-nd ^ — ^7,0 
are inverses of each other. 

Theorem 2.5. ^(<I>(/)) = / for any f € MVg and <I>(^(F)) = F for any F € J^g- 
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Proof. First, if / G MVg F = <!>(/), then by Prop. 12. li the roots of / are exactly the minimal 
vertices in each tree component of F. Those in turn are roots for ^{F). In applying algorithm 
B to F, we note that the order vr = viV2 ■ ■ - Vn is exactly the order in which vertices of G will be 
processed when running algorithm A on /. That is, Pi = {vi, . . . ,Vi}, and fj+i is not a root of /, 
then Qi is the set of vertices which are adjacent (via edges in F) to those in Pi. By the construction 
of algorithm A, a vertex w is found by v if and only if there are f{w) many edges connecting w to 
vertices that are processed before v, or equivalently, to vertices u with tt~^{u) < ■k~^{v). Since in 
$(/), V = wP, we have ^'(^>(/)) = /. 

Conversely, we prove that ^{^{F)) = F hy showing that $(^'(-F)) and F have the same set of 
edges. First note that the minimal vertices of the tree components of F are exactly the roots of 
/ = ^{F), which then are the minimal elements of trees in ^{f). Edges of F are of the form {v, v^}, 
where v is not a minimal vertex in its tree component. We now show that when applying algorithm 
A to ^(F), vertex v is found by v^. Note that f{v) = \{vj\{v,Vj) G E{G),tt~'^ (vj) < tt~'^{vP)}\. In 
the implementation of algorithm A, the valuation on v drops by 1 for each adjacent vertex that is 
processed before v. When it is w^'s turn to be processed, vali{v) drops from to —1. Thus finds 
V, and {v, v^} is an edge of 

□ 

Since the roots of the G-multiparking function correspond exactly to the minimal vertices in 
the tree components of the corresponding forest, in the following we will refer to those vertices as 
roots of the forest. 

3 Examples of the bijections 

The bijections ^7,g and ^7,G' as defined above via algorithms A and B, allow a good deal of 
freedom in implementation. In algorithm A, as long as 7 is well-defined at every iteration of Step 
2, one can obtain fa^j+i, Pj+i, Qi+i and Fj+i and proceed. Recall that 7 is a function from W, the 
set of ordered pairs (F, VF), to V{G) such that 7(F, W) € W , where F is any sub-forest of G (not 
necessarily spanning) and 1^ is a non-empty subset of Leaf{F) or consists of an isolated point of 
F. 

When restricting to G-parking functions, (i.e., G-multiparking functions with only one root), 
the descriptions of the bijections <I> and ^ are basically the same as the ones given by Chebikin and 
Pylyavskyy |2j, where the corresponding sub-structures in G are spanning trees. However our family 
of bijections, each defined on a choice function 7, is more general than the ones in [2], which rely on 
a proper set of tree orders. A proper set of tree orders is a set n(G) = {vr(T) : T is a subtree of G} 
of linear orders on the vertices of T, such that for any -y G T, w <-k{t) "^^i and if T' is a subtree 
of T containing the least vertex, tt{T') is a suborder of tt(T). Our algorithms do not require there 
to be a linear order on the vertices of each subtree. In fact, for a spanning tree T of a connected 



8 



graph G, the proper tree order tt{T), if it exists, must be the same as the one defined in Step 1 
of algorithm B. But in general, for two spanning trees T and T' with a common subtree t, the 
restrictions of vr(T) and vr(T') to vertices of t may not agree. Hence in general the choice function 
cannot be described in terms of proper sets of tree orders. In addition, our description of the map 

in terms of a dynamic process, provides a much clearer way to understand the bijection, and 
leads to a natural classification of the edges of G which plays an important role in connection with 
the Tutte polynomial (c.f.§4). 

Different choice functions 7 will induce different bijections between A4Vg and Tq- In this sec- 
tion we give several examples of choice functions that have combinatorial significance. In Example 
1 we explain how to translate a proper set of tree orders into a choice function. Hence the family 
of bijections defined in [2j can be viewed as a subfamily of our bijections restricted to G-parking 
functions. The next three examples have appeared in We list them here for their combinato- 
rial significance. Example 5 is the combination of breadth-first search with the Q-sets equipped 
with certain data structures. It is the one used to establish connections with Tutte polynomial 
in §4. The last example illustrates a case where 7 cannot be expressed as a proper set of tree 
orders. We illustrate the corresponding map ^-y^c for examples 2-6 on the graph G in Figure Q A 
G-multiparking function / is indicated by ^^i/ f{i)" on vertices, where i is the vertex label. 

1/(30 3/0 




7/0 

Figure 1: A graph and a multiparking function. 

In each example, we will show the resulting spanning forest by darkened edges in G. Again 
each vertex will be labeled by a pair i/j, where i is the vertex labels, and j = valn{i), where n = 7. 
Beneath that, a table will record the sets Qt and Pt for each time t. In each Qt, the vertex listed 
first is the next to be processed. 
Example 1. 7 with a proper set of tree orders. 

We define the choice function that corresponds to a proper set of tree orders. Here we should 
generalize to the proper set of forest orders, i.e., a set of orders vr(F), defined on the set of vertices 
for each subforest F of G, such that for any v £ F, v <n{F) and if F' is a subforest of F with 
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the same minimal vertex in each tree component, t^{F') is a suborder of tt{F). In this case, define 
7(-F, W) = V where v is the minimal element in W under the order 7r{F). Examples 2-4 are special 
cases of this kind. 

Example 2. 7 with a given vertex ranking. 

Given a vertex ranking a Sn define "YaiFjW) := v, where v is the vertex in W with minimal 
ranking. In particular, if a is the identity permutation, then the vertex processing order is the 
vertex- adding order of 0. In this case, in Step 2 of algorithm A, we choose v to be the least vertex 
in Qi-i and process it at time i. The output of algorithm A is 



1/00 3/-2 




7/-2 



The Qi and Pi for this instance are as follows. 



t 





1 


2 


3 


4 


5 


6 


7 


Qt 


{1} 


{2,3} 


{3,6} 


{4,6} 


{5,6} 


{6,7} 


{7} 





Pt 





{1} 


{1,2} 


{1,2,3} 


{1,2,3,4} 


{1,2,3,4,5} 


{1,2,3,4,5,6} 


{1,2,3,4,5,6,7} 



Example 3. 7 with depth-first search order. 



The depth-first search order is the order in which vertices of a forest are visited when performing 
the depth-first search, which is also known as the preorder traversal. Given a forest F with tree 
components Ti, T2, . . . , T^, where 1 = ri < r2 < • • • < are the corresponding roots, the order <df 
is defined as follows. (1) For any v & Ti, w £ Tj and i < j, v <df w. (2) For any v 7^ r^, Kdf v. 
(3) If = and v < w, v <df w. (4) For any v, let F[v] be the subtree of F rooted at v. If 
V € w € F[w'] and v' <fif w' , then v <df w. For example, the depth-first search order on the 

below tree is 1 <df 2 <df 3 <df 6 <df 4 <df 5. 




Figure 2: A tree with 6 vertices. 
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The choice function ^df with depth- first search order is then defined as "ydf{F, W) = v where v 
is the minimal element of W under the depth-first search order <df of F. Here is the output of 
algorithm A with the choice function ^^f on the example in Figure ^ 

l/oo 3/-5 




7/-2 

The Qi and Pj for this instance are as follows. 



t 





1 


2 


3 


4 


5 


6 


7 


Qt 


{1} 


{2,3} 


{6,3} 


{4,3,7} 


{5,3,7} 


{7,3} 


{3} 





Pt 





{1} 


{1,2} 


{1,2,6} 


{1,2,4,6} 


{1,2,4,5,6} 


{1,2,4,5,6,7} 


{1,2,3,4,5,6,7} 



Example 4. 7 with breadth-first search order. 

Breadth-first search is another commonly used tree traversal in computer science. Given a forest 
F, whose tree components are Tj with roots rj, (1 < i < k), and 1 = ri < r2 < • • • < r^, the order 
<bf is defined as follows. (1) For any v & Ti, w £ Tj and i < j, v <{,/ w. (2) Within tree Tj, for 
each V Ti, let height hTi{v) of v be the number of edges in the unique path from v to the root 
Tj. We set V <bf w if hxiiv) < /^T^(^«), or else if hTi{v) = hxiiw) and v < w. For example, the the 
breadth-first search order for the tree in Figure [21 is 1 <bf 2 <hf 4 <f,y 3 <5j 5 <bf 6. 

The choice function 7^,^ with breadth-first search order is defined as 7b/ (-F, W) = v where v is 
the minimal element of W under the breadth-first search order <bf of F. Here is the output of 
algorithm A with the choice function -f{,f on the example in Figure ^ 



1/00 3/-2 




7/-2 
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The Qi and Pi for this instance are as follows. 



t 





1 


2 


3 


4 


5 


6 


7 


Qt 


{1} 


{2,3} 


{3,6} 


{4,6} 


{6,5} 


{5,7} 


{7} 





Pt 





{1} 


{1,2} 


{1,2,3} 


{1,2,3,4} 


{1,2,3,4,6} 


{1,2,3,4,5,6} 


{1,2,3,4,5,6,7} 



Example 5. Breadth- first search with a data structure on Qj. 

In this case, new vertices enter the set Qi in a certain order, and some intrinsic data structure 
on Qi decides which vertex of Qi is to be processed in the next step. A typical example is that 
of breadth-first search with a queue, in which case each Qi is an ordered set, (i.e., the stage of a 
queue at time i). New vertices enter Qi in numerical order, and 7 chooses the vertex that entered 
the queue earliest. 

This example can also be defined by a modified breadth-first search order, which we call breadth- 
first order with a queue, and denote by <bf,q- Given a forest F, whose tree components are Tj with 
root Tj, (1 < i < k), and 1 = ri < r2 < • • • < r^, the order <bf,g is defined as follows. (1) For any 

V G Ti, w & Tj and i < j, v <bf,q w. (2) Within tree Tj, the root is minimal under <bf,q- (3) 

V <bf,q vj ii yP <bf,q w^. (4) If = and y < y <bf,q w. For example, the breadth-first search 
order with a queue for the tree in Figured is 1 <bf,q 2 <bf,q 4 <bf,q 3 <bf,q 6 <bf,q 5. 

The choice function 7 associated with this order is denoted by 7b/.q, and is used in §4. The 
following is the output of algorithm A with 7fe/,g on the graph in Figure ^ 



1/00 3/-2 




7/-2 

The Qi and Pi for this instance are as follows, where each Qi is an ordered set, and the first element 
in Qi is the next one to be processed. 



t 





1 


2 


3 


4 


5 


6 


7 


Qt 


(1) 


(2,3) 


(3,6) 


(6,4) 


(4,5,7) 


(5,7) 


(7) 





Pt 





{1} 


{1,2} 


{1,2,3} 


{1,2,3,6} 


{1,2,3,4,6} 


{1,2,3,4,5,6} 


{1,2,3,4,5,6,7} 



Another typical structure is to let Qi be the stage of a stack at time i, that is, it pops out the 
vertex that last entered. We can also combine the other vertex orders with a queue or stack for the 
Q-sets. 
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Example 6. A choice function 7 that cannot be defined by a proper set of tree orders. 
Let 




the second minimal vertex of W, if \W\ > 2. 




Then the order on the left tree is 156342, and the one on the right tree is 153462, which do not 
agree on the subtree consisting of vertices 1356. Hence it can not be defined via a proper set of 
tree orders. 



Remark. Note that the bijection given in [2| applies to a general directed graph. We explain 
how our algorithms could be slightly modified to apply to that case, too. Let D denote a general 
directed graph on [n]. An oriented spanning forest F of D is a subgraph of G such that (1) the 
edges, when ignoring the orientation, do not form a cycle, and (2) for each v in a tree component 
with minimal vertex r, there is a unique directed path from -y to r. Again we denote by v'^ the 
vertex lying on the directed path from f to r with {v, v^) G E{D). We say that the minimal vertex 
in each tree component of F is a root of F. 

The definition of a D-multiparking function is the same as that of the undirected multiparking 
function, except that outdegu{i) is the number of edges going from i to vertices not in U. Again 
we say the vertices v with g{v) = 00 are the roots of the Z)-multiparking function g. 

Spanning forests do not contain loops, as loops are a trivial kind of cycle. Hence we can assume 
D is loop less without loss of generality. We allow D to have multiple edges. But to distinguish 
between multiple edges of -D, we fix a total order on the set of edges going from i to j, for each 



The maps <^ and ^ can be modified accordingly to give a bijection between the set of D- 
multiparking functions to the set of oriented spanning forests of D, which carry the roots of mul- 
tiparking functions to the roots of spanning forests. The only modifications we need to make 
are: 

For Algorithm A. 

In Step 3, lower the value of u; by 1 for each directed edge from tt; to if it; is not a root. Load 
w to Qi whenever vali{w) < 0. For each such add the {k + l)-st edge between from w to v ii 
vali-i{w) = k > 0. 
For Algorithm B. 

In Step 1. Let W be the set {v ^Vi -.^w e Vi such that {v, w) G E{D)}. 




i + J- 



13 



In Step 2. For any nonroot vertex v lying in the tree with root r^,, and v ^ r„, set f{v) to be 
k + Vj) e E{G), TT~^{vi) < 'K~^{v^)}\ if the edge {v, v^) in F is the (k + l)-st edge in the set 

of edges from v to in D. 

These modified algorithms for directed graphs cover the case for an undirected graph G, provided 
that one views G as a digraph, where each edge {u, v} of G is replaced with two directed edges 
{u, v) and {v, u). 

The notion of G-parking functions, as proposed in ^Sj, is closely related to the critical configu- 
rations of the chip-firing games, (also known as sand-pile models). The generalization of chip-firing 
games with multiple sources is given in [H], where they are called Dirichlet games. In the first 
author shows that modified G-multiparking functions (which in Condition (A), instead of requiring 
that i = min{U), one requires i to belong to a prefixed subset of vertices), are the corresponding 
counterpart for critical configurations of the Dirichlet games. In fact, both are in one-to-one corre- 
spondence with the set of rooted spanning forests of G, as well as a set of objects called descending 
multitr aver sals on G. 

4 External activity and the Tutte polynomial 

4.1 F-redundant edges 

A forest F on [n] may appear as a subgraph of different graphs, and a vertex function / may be 
a G-multiparking function for different graphs. In this section we characterize the set of graphs 
which share the same pair {F, /). Again let G be a simple graph on [n], and fix a choice function 
7. For a spanning forest F of G, let / = ^'-y^G'(F)- We say an edge e of G — F is F -redundant if 
^7,G-{e}(-^) = /• Note that we only need to use the value of 7 on {F\ W) where F' is a sub-forest 
of F. Hence G'_{e}(F) is well-defined. 

Let vr be the order defined in Step 1 of Algorithm B. Note that vr only depends on F, not the 
underlying graph G. Recall that denotes the parent vertex of vertex v in some spanning forest. 
We have the following proposition. 

Proposition 4.1. 

An edge e = of G is F -redundant if and only if e is one of the following types: 

1. Both V and w are roots of F. 

2. V is a root and w is a non-root of F, and 7r~^{w) < 7r~^[v). 

3. V and w are non-roots and 7r~^(t;^) < 7r~"'^('u;) < 7r~^(t)). In this case v and w must lie in the 
same tree of F. 
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Proof. We first show that each edge of the above three types are F-redundant. Since for any root 
r of the forest F, f{r) = oo, the edges of the first two types play no role in defining the function /. 
And clearly those edges are not in F. Hence they are F- redundant. 

For edge {v.,w) of type 3, clearly it cannot be an edge of F. Since f{v) = 4j^{vj\{v.,Vj) G 
E{G),'K~^ (vj) < 7r~^(uP)}, and 7r~^{w) > tt~^{vP), removing the edge {v,w} would not change the 
value of f{v). This edges has no contribution in defining f{u) for any other vertex u. Hence it is 
F-redundant. 

For the converse, suppose that e = {v, w} is not one of the three type. Assume w is processed 
before v in vr. Then v is not a root, and w appears before v^. Then removing the edge e will change 
the value of f{v). Hence it is not F-redundant. 

□ 

Let Ri{G; F), R2{G; F), and i?3(G; F) denote the sets of F-redundant edges of types 1, 2, and 
3, respectively. Among them, i?3(G;F) is the most interesting one, as Ri{G; F) and R2{G;F) are 
a consequence of the requirement that /(r) = oo for any root r. Let R{G; F) be the union of these 
three sets. Clearly the F-redundant edges are mutually independent, and can be removed one by 
one without changing the corresponding G-multiparking function. Hence 

Theorem 4.2. Let H be a subgraph of G with V{H) = V{G). Then ^^,g{F) = ^'y,H{F) if and 
only ifG- R{G; F) CH CG. 

4.2 A classification of the edges of G 

The notion of F- redundancy allows us to classify the edges of a graph in terms of the algorithm A. 
Roughly speaking, the edges of any graph can be thought of as either lowering val{v) for some v 
to 0, being in the forest, or being F-redundant. Explicitly, we have 

Proposition 4.3. 

Let f be a G-multiparking function and let F = •&(/). Then 



\E{G)\ = 




Proof. For each non-root vertex v, the number of different values that vali{v) takes on during the 
execution of algorithm A is f{v) + 1 + ny, where n„ = —valn{v). At the beginning, valQ{v) = f{v). 
The value vali{v) then is lowered by one whenever there is a vertex w which is adjacent to v and 
processed before v'p. When is being processed, vali{v) = —1, and the edge {v^, v} contributes to 
the forest F. Afterward, the value of vali{v) decreases by 1 for each F-redundant edge {n, f} with 
'K~^{u) < 'K~^{v). Summing over all non-root vertices gives 

E '^^5<.w = E /(^) + 1^(^)1+ E 

v:f{v)y^oo v:f{v)j^oo v:f{v)^oo 
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where deg<^{v) = \{{w,v} G E{G)\7r~'^ (w) < tt~^{v)}\. 

The edges that lower val{v) below —1 are exactly the F-redundant edges of type (3) in Prop. 
EH hence E«:/(«)^oo = 1^3 (G; F)\. On the other hand, Y.v:f{v)^oo deg<^{v) is exactly \E{G) \ - 
\Ri{G;F)\ - |i?2(G;F)j. The claim follows from the fact that the sets Ri{G; F), R2{G; F), and 
Rs{G;F) are mutually exclusive. 

□ 

One notes that for roots of / and F = $(/), \Ri{G; F)\ + \R2{G; F)\ is exactly Y^root v ^e5<^('u), 
where tt is the processing order in algorithm A. But it is not necessary to run the full algorithm 
A to compute |iii(G;-F)| + |i?2(G;-F)|. Instead, we can apply the burning algorithm in a greedy 
way to find an ordering tt' = v[v'2 ■ ■ ■ v'^ on V{G): Let v'l = 1. After determining v'^^, . . . ,v[_i, if 
Vi = V{G) — {ui, . . . ,f •__]^} has a well-behaved vertex, let v[ be one of them; otherwise, let v[ be 
the minimal vertex of 1^, (which has to be a root.) 

vr' may not be the same as vr, but they have the following properties: 

1. Let ri < r2 < • • • < Tfc be the roots of /. Then ri, r2, . . . , appear in the same positions in 
both vr and vr'. 

2. The set of vertices lying between rj and rj+i are the same in vr and vr'. In fact, they are the 
vertices of the tree Tj with root rj in F = <!*(/). 

It follows that for any root vertex v, deg<^^{v) = deg<^^,{v). The value of deg^^{v) {v root) 
can be characterized by a global description: Let lA^ be the collection of subsets U of V{G) such 
that V = min(?7), and U does not have a well-behaved vertex. lAv is nonempty for a root v since 
U = {v} \s such a set. Then 

deg^^{v) = min outdegu{v). 

U G Vlv 

We call deg^^{v) the record of the root -y, and denote it by rec{v). Then 
\Ri{G;F)\ + \R2{G;F)\= ^ deg<^,{v) = rec{v) 

root V root v 

is the total root records. Let Rec{f) = |-Ri(G; F)\ + \ R2{G; F)\. It is the number of F-redundant 
edges adjacent to a root. By the above greedy burning algorithm, the total root records Rec{f) 
can be computed in linear time. 

4.3 A new expression for Tutte polynomial 

In this subsection we relate G-multiparking functions to the Tutte polynomial tG{x,y) of G. We 
follow the presentation of [H] for the definition of Tutte polynomial and its basic properties. Al- 
though the theory works for general graphs with multiedges, we assume G is a simple connected 
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graph to simplify the discussion. There is no loss of generality by assuming connectedness, since 
for a disconnected graph, tcix^y) is just the product of the Tutte polynomials of the components 
of G. We restrict ourselves to connected graphs to avoid any possible confusion when we consider 
their spanning forests. The modification when G has multiple edges is explained at the end of §3. 

Suppose we are given G and a total ordering of its edges. Consider a spanning tree T of G. An 
edge e S G — T is externally active if it is the largest edge in the unique cycle contained in T U e. 
We let 

£A{T) = set of externally active edges in T 

and ea{T) = \£A{T)\. An edge e € T is internally active if it is the largest edge in the unique 
cocycle contained in (G — T) U e. We let 

ZA{T) = set of internally active edges in T 

and ia{T) = \IA{T)\. Tutte |2J then defined his polynomial as 

iG(x,y)= (1) 

TCG 

where the sum is over all spanning trees T of G. Tutte showed that tc is well-defined, i.e., inde- 
pendent of the total ordering of the edges of G. Henceforth, we will not assume that the edges of 
G are ordered. 

Let H he a (spanning) subgraph of G. Denote by c{H) the number of components of H. Define 
two invariants associated with H as 

aiH) = c{H) - 1, a\H) = \E{H)\ - \V{G)\ + c{H). (2) 
The following identity is well-known, for example, see PP. 
Theorem 4.4. 

tG(l + x,l + y)= (3) 

H<ZG 

where the sum is over all spanning subgraphs H of G. 

Recall that the breadth-first search (BFS) is an algorithm that gives a spanning forest in the 
graph H. Assume V{G) = [n]. We will use our favorite description to express the BFS as a queue 
Q that starts at the least vertex 1. This description was first introduced in to develop an exact 
formula for the number of labeled connected graphs on [n] with a fixed number of edges, and was 
used by the second author in to reveal the connection between the classical parking functions 
(resp. /c-parking functions) and the complete graph (resp. multicolored graphs). 
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Given a subgraph H of G with V{H) = V{G) = [n], we construct a queue Q. At time 0, Q 
contains only the vertex 1. At each stage we take the vertex x at the head of the queue, remove x 
from the queue, and add all unvisited neighbors ui, . . . ,ut^ of x to the queue, in numerical order. 
We will call this operation "processing x" . If the queue becomes empty, add the least unvisited 
vertex to Q. The output F is the forest whose edge set consists of all edges of the form {x,Ui} 
for i = 1, . . .tx- We will denote this output as -F = BFS{H). Figure |H1 shows the spanning forest 
found by BFS for a graph G. 




Figure 3: Spanning forest found by BFS. 



The queue Q for Figure IHl is 



t 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Q 


(1) 


(3,4) 


(4,8) 


(8,7) 


(7) 


(6,9) 


(9) 


(2) 


(5,10) 


(10,11) 


(11) 






For a spanning forest F of G, let us say that an edge e G G — F is BFS- externally active 
if BFS{F D e) = F. A crucial observation is made by Spencer ^U]: An edge {v,w} can be 
added to F without changing the spanning forest under the BFS if and only if the two vertices 
V and w have been present in the queue at the same time. In our example of Figure |SJ edges 
{3, 4}, {4, 8}, {7, 8}, {6, 9}, {5, 10} and {10, 11} could be added back to F. We write S{F) for the 
set of BFS-externally active edges. 

Proposition 4.5 (Spencer). If H is any subgraph and F is any spanning forest of G then 

BFS{H) = F if and only if F C H C F U £{F). 

Now consider the Tutte polynomial. Note that if BFS{H) = F, then c{H) = c{F). So 
a{H) = c{F) - 1 and a*{H) = \E{H) \ - \E{F)\ = \£{F) n H\. Hence if we fix a forest F and sum 
over the corresponding interval [F, F U £{F)]^ we have 

^ ^'^{H)y.'iH) ^^ciF)-l yl^l=x<^(^)-l(l + y)l^(^)l 

H:BFS(H)=F ACe{F) 
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Summing over all forests F, we get 

HCG FCG 

Or, equivalently, 

tG{l + x,y)=J2x<^^-'y\^(^K (4) 
FCG 

To evaluate £{F), note that when applying BFS to a graph H, the queue Q only depends on 
the spanning forest F = BFS{H). Given a forest F, the processing order in Q is a total order 
<Q=<Q {F) on the vertices of F satisfying the following condition: Let Ti,T2, ■ ■ ■ ^T^ be the tree 
components of F with minimal elements ri = 1 < r2 < • • • < r^. Then (1) If is a vertex in tree Tj, 
w is a vertex in tree Tj and i < j, then v <q w. (2) Among vertices of each tree Tj, is minimal 
in the order <q. (3) For two non-root vertices v,w in the same tree, v <q w if <q w^. In the 
case = w^, v <q w whenever v < w. 

Comparing with the examples in §3, we note that <q is exactly the order <bf,q described in 
Example 5 of §3, as breadth-first order with a queue. Fix the choice function 7 = 76/,g, the one 
associated to <bf,q and consider the maps ^*7,g and ^I'7,G• Given F, the condition that two vertices 
V, w have been present at the queue Q at the same time when applying BFS to F is equivalent to 

<bf,q w <bf,q V or <bf,q V <bf,q w. That is, an edge is BFS-externally active if and only if it 
is an F-redundant edge of type 3, as defined in §4.1. It follows that £{F) = Rz{G; F). 

Therefore by Prop. 14.31 

\£{F)\ = \R,{G;F)\ = \E{G)\-\E{F)\-( f{v))-Rec{f), 

where / = ^^^g{F^ is the corresponding G-multiparking function. Note that \E(F')\ = n — c{F), 
and c{F) = r{f), where r(/) is the number of roots of /. Therefore 

Theorem 4.6. 

toil + x,y) = J2 ^(^)) , 

/ 

where the sum is over all G-multiparking functions. 

For a G-multiparking function /, where G is a graph on n vertices, we call the statistics 
|£^(G)| — n + r{f) — Rec{f) — Ylvf(v)^oo f(^) reversed sum of /, denote by rsum{f). The 
name comes from the corresponding notation for classical parking functions, see, for example, jl5j . 
Theorem 14.61 expresses Tutte polynomial in terms of generating functions of r(/) and rsum[f). In 
[U] Gessel and Sagan gave a similar expression, in terms of £dfs{F)^ the set of greatest-neighbor 
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externally active edges of F, which is defined by applying the greatest-neighbor depth-first search 
on subgraphs of G. Combining the result of [HI (Formula 5), we have 

XtGil + X,y) = = ^ = ^ :,rif)yrsumif)_ (5) 

FCG FCG feMVa 

That is, the three pairs of statistics, (c(F), \£j:)FsiF)\) and (c(F), for spanning forests, 

and {r{f),rsum{f)) for G-multiparking functions, are equally distributed. 

Remark. Alternatively, one can prove Theorem l4.6l bv conducting neighbors-first search (NFS), 
a tree traversal defined in |UJ §6], and using 7 = ^df, the choice function associated with the depth- 
first search order, (c.f. Example 3, §3). Here the NFS is another algorithm that builds a spanning 
forest F given an input graph H. The following description is taken from [5]. 

NFSl Let F = 0. 

NFS2 Let v be the least unmarked vertex in V and mark v. 

NFS3 Search v by marking all neighbors of v that have not been marked and adding to F all edges 
from V to these vertices. 

NFS4 Recursively search all the vertices marked in NFS3 in increasing order, stopping when every 
vertex that has been marked has also been searched. 

NFS5 If there are unmarked vertices, then return to NFS2. Otherwise, stop. 

The NFS searches vertices of in a depth-first manner but marks children in a locally breadth- 
first manner. Figure 0] shows the result of NFS, when applies to the graph on the left of Figure 

121 




6 9 

Figure 4: Spanning forest found by NFS. 



Similarly, one defines SisfFsiF")^ the set of edges externally active with respect to NFS, to be 
those edges e £ G - F such that NFS{F Ue) = F. Then Prop. EH and Eq. ^ hold again when 
we replace BFS with NFS, and £{F) with £nfs{F). 
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Now let 7 = 7rf/ and use the bijections 07^y,G and ^^^^.^g, one notices again that an edge is 
externally active with respect to NFS if and only if it is -F-redundant of type 3. And hence we get 
another proof of Theorem 14.61 

An interesting specialization of Theorem 14.61 is to consider tG{l,y), the restriction to spanning 
trees of G and G-parking functions. For a G-parking function /, or equivalently a G-multiparking 
function with exactly one root (which is vertex 1), r(/) = 1 and Rec{f) = 0. Hence rsum{f) = 
\E{G)\ - n + 1 - f{v). Thus we obtain 

f : G-parking functions 

An equivalent form of this result, in the language of sand-pile models, was first proved by Lopez 
jl6j using a recursive characterization of Tutte polynomial. A bijective proof was given by Cori and 
Le Borgne in |3] by constructing a one-to-one correspondence between trees with external activity 
i (in Tutte's sense) to recurrent configurations of level i, which is equivalent to G-parking functions 
with reversed sum i. Our treatment here provides a new bijective proof. 

In |n] it is shown that, restricted to simple graphs, the greatest-neighbor externally active edges 
of F are in one-to-one correspondence with certain inversions of F. For a simple graph G, view 
each tree T of F as rooted at its smallest vertex. An edge {u, v} is greatest-neighbor externally 
active if and only if v is a descendant of u, and w > v where w is the child of u on the unique u — v 
path in F, (that is, u = w^). Call such a pair {w,v} a G-inversion. And denoted by Ginv{F) the 
number of G-inversions of the forest F. Then we have the following corollary. 

Corollary 4.7. Let Tk{G) he the set of spanning forests of G with exactly k tree components. And 
MVk{G) he the set of G-multiparking functions with k roots. Then 



V ^ Ginv{F) ^ ^ ,,i'''sum{f) 



^ y ' ' = 2^ y 



In particular, when G is the complete graph -ftTn+i and /c = 1, we have the well-known result 
on the equal-distribution of inversions over labeled trees, and the reversed sum over all classical 
parking functions of length n, (for example, see |13[ l2Uj ) 

T on [n+1] QSPn 

where Pn is the set of all (classical) parking functions of length n. 



5 Enumeration of G-multiparking functions and graphs 

In this section we discuss some enumerative results on G-multiparking functions and substructures 
of graphs. 
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Theorem 5.1. The number of G-multiparking functions with k roots equals the number of spanning 
forests of G with k components. In particular, for connected graph G, the number of G-multiparking 
functions is Tq{2, 1). Among them, those with an odd number of roots is counted by ^[Tg{2, 1) + 
Tq{0, 1)), and those with an even number of roots is counted by ^(Tg'(2, 1) — Tg{0, 1)). 

Proof. The first two sentences follow dhectly from the bijections constructed in §2, and Theorem 
14.61 For the third sentence, just note that Tg(0, 1) = ^f{—^Y^^^~^ is the difference between the 
number of G-multiparking functions with an odd number of roots, and those with an even number 
of roots. □ 

Another consequence of Theorem 14.61 and its proof is an expression for the number of spanning 
subgraphs with a fixed number of components and fixed number of edges, in terms of (BFS)- 
external activity and G-multiparking functions. It is a generalization of the expectation formula in 
[T^ . which is the special case for complete graph Kn- 

Theorem 5.2. Let G be a connected graph. The number ^t,k{G) of spanning subgraphs H with t 
components and V{G) — 1 + k edges is given by 

where the first sum is over all spanning forests with t components, and the second sum is over all 
G-multiparking functions with t roots. 

Proof. For any spanning forest F with k components, the number of spanning subgraphs H with 
V{G) -1 + k edges such that BFS{H) = F is given by {^^p) . □ 

Next we give a new expression of the tK„_f_i{x,y) in terms of classical parking functions. It 
enumerates the classical parking functions by the number of critical left-to-right maxima. Given a 
classical parking function b = . . . , bn), we say that a term bi = j is critical if in b there are 
exactly j terms less than j, and exactly n — l—j terms larger than j. For example, in b = (3, 0, 0, 2), 
the terms hi = 3 and 64 = 2 are critical. Among them, only 61 = 3 is also a left-to-right maximum. 

Let a(b) be the number of critical left-to-right maxima in a classical parking function b. We 
have 

Theorem 5.3. 

hePn 

where Pn is the set of classical parking functions of length n. 

Proof. Let F be a spanning forest on [?7- + 1] with tree components Ti, . . . , T^, where Tj has minimal 
vertex rj, and ri < r2 < ■ ■ ■ < r^. We define an operation merge{F) which combines the trees 
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Ti, . . . ,Tk by adding an edge between rj with Wi-i for each i = 2, k, where Wi-i is the vertex 
of Tr-i that is maximal under the order <bf,q- Denote by Tp = merge{F) the resulting tree. We 
observe that for the forest F and the tree Tp, the queue obtained by applying BFS are exactly the 
same. This implies that F and Tp have the same set of BFS-externally active edges. 

Conversely, given T and an edge e = {w^v} € T where w <bf,q v. We say the edge e is critical 
in T if merge{T \ {e}) = T. Assume T \ {e} = Ti U T2 where w £ Ti and v € T2. By the definition 
of the merge operation, e is critical if and only if w is the maximal in Ti under the order <bf,q, and 
V is vertex of the lowest index in T2. In terms of the queue obtained by applying BFS to T, it is 
equivalent to the following two conditions: (1) There is a set Qi such that Qi = {v}, and v does 
not belong to any other Qi. (2) v is of minimal index among the set of vertices processed after v. 

Consider the maps $7,0 and ^'7,0 with 7 = 'ybf,q and G = Kn+i- Let / = ^'7^g'(T), and write / 
as a sequence (/(2), /(3), . . . , /(n+1)). (There is no need to record /(I), as /(I) = 00 always.) Then 
an edge {w, v} is critical in T if and only if (1) f[v) is critical in the sequence (/(2), . . . , /(n + 1)), 
and [2) w > V for any vertex w with f{w) > fiv). That is, f{v) is a left-to-right maximum in the 
sequence (/(2), /(3), . . . , /(n + 1)). 

Now fix a spanning tree T of Kn+i and let Merge{T) be the set of spanning forests F such 
that merge{F) = T. Then an G Merge{T) can be obtained from T by removing any subset A 
of critical edges, in which case c{F) = c{T) + \ A\ = 1 -|- \ A\. This, combined with the fact that 
£{F) = £{T), gives us 

F&Merge{T) A 

where A ranges over all subsets of critical edges of T. Under the correspondence T ^ f = ^^^ciT) 
and considering / as a sequence (/(2), . . . -|- 1)), j<f(T)| is just (2) — Y17='2 critical 
edges of T correspond to critical left-to-right maxima of the sequence. Hence the sum in © equals 

yi^mi(i + ^)"{/T) = (i + ^)«(/T)y©-EriV». 

Theorem 15.31 follows by summing over all trees on [n -1- 1] . □ 

Finally, we use the breadth-first search to re-derive the formula for the number of subdigraphs 
of G, which was first proved in jH] using DFS, and extend the method to derive a formula for the 
number of subtraffics of G. 

Let G be a graph. A directed subgraph or subdigraph of G is a digraph D that contains up to 
one copy of each orientation of every edge of G. Here for an edge {n, v} of G we permit both (u, v) 
and {v,u) to appear in a subdigraph. 

For any subdigraph D of G, we apply the BFS to get a spanning forest of D. The only difference 
from the subgraph case is that when processing a vertex x, we only add those unvisited vertices u 
such {x, u) is an edge of D. 
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If digraph D has BFS forest F, write -^+(1?) = F. Note that we can view F as an oriented 
spanning forest, where each edge is pointing away from the root (i.e., the minimal vertex) of the 
underlying tree component. Say a directed edge e ^ F is directed BFS externally active with respect 
to F if T'^{F L) e) = F. Denote by £'^{F) the set of directed BFS-externally active edges. Then 
we have the following basic proposition, which is the analog in the undirected case. 

Proposition 5.4. // D is any subdigraph and F is any spanning forest of G then 

f+{D) = D if and only if F C D Q F U £'+(F). 

Now we characterize the directed BFS-externally active edges by the set £{F), the BFS- 
externally active edges for the undirected graph G. Let {u,v} be an edge of G with u <bf,q v. 
If {u,v} G E{F), then the backward edge {v,u) can be added without changing the result of (di- 
rected) breadth-first search, that is, {v,u) € £~^{F). If {u,v} € £{F), then both {u,v) and {v,u) 
are in £~^{F). If {u,v} is not in the forest F or £{F), then {v,u) is in £^{F). Together we have 

\E{G)\ = \£+{F)\-\£{F)\. 

Therefore 

Theorem 5.5. If G has n vertices, then 

Y: = xy--\l + y)l^(^)l tail + -, 1 + y), (7) 

where the sum is over all subdigraphs of G. 
Proof. 

^ ^ D:f+{D)=F 

= 5]x^(^V^(^)l(l+y)l^^(^)l 
F 

= ?/"(i+2/)i^(^)iE(-) 

p \y / 

= xy"-i(l + y)l^(«)l tG{l + -A + y). 

V 

□ 

Next we consider a slightly complicated problem. The sub-traffic K of G, where K \s a, partially 
directed graph on 1^(G), is obtained from G by replacing each edge {u,v} of G by (a) 0, (b) a 
directed edge {u,v), (c) a directed edge {v,u), (d) two directed edges {u,v) and {v,u), or (e) an 
undirected edge {u, v}. We proceed as we did before. For each subtraffic K, we apply the directed 
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breadth- first search to get a spanning forest F: The queue starts with the minimal vertex 1. At 
each iteration, we take the vertex x at the head of the queue, remove x from the queue, and add all 
unvisited vertices u if {x,u) S E{K) or {x,u} G E{K). Add the directed edge {x,u) to the forest 
F if {x,u) € E{K). Otherwise, add the undirected edge {x,n} to F. The output is a forest [n] in 
which each edge is either a directed edge oriented away from the minimal vertex of the underlying 
tree, or an undirected edge. Let A be the set of directed edges. Denote by {F, A) the output forest 
and write BFS{K) = (F, A). Note that (F, A) is itself a sub-traffic of G. 

Given a pair {F,A) with directed edges A C E{F), we have the following characterization of 
edges that can be added to {F,A), without changing the BFS result, (i.e., BFS{{F, A) U e) = 
BFS{F,A).) 

1. For each directed edge {u, v) in A Q E{F), we can add back {v, u) without changing the result 
of the spanning forest. 

2. For each BF-externally active edge {u,v} of F, we can add back any one of {u,v), {v,u) and 
{u,v}, or both {u,v) and {v,u) at the same time. 

3. For each edge not in F U £{F), we can add back one of the undirected edge {u,v} and the 
direct {u, v) if u is processed after v in the queue. 

There is no further restriction on how the edges can be added back in addition to the above 
mentioned cases. Then we have 

Theorem 5.6. Let G be a connected graph. Then 

where the sum is over all subtraffic of G. 
Proof. 

K F ACE{F) K:BFS{K)={F,A) 

where F is over all spanning forests of G, and A is a subset of the edges of F. A subtraffic K has 
BFS{K) = (F, A) if and only if it is obtained from F by adjoining some edges as described in the 
preceding three cases. Considering the contribution of each type, we have 

K:BFS{K)={F,A) 
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Hence 

/ \ \E{F)\ /I , Q , 2\ \£iF)\ 

K F \ ^ y / \-ry / ACE{F) 

/ „, \\E(F)\ ,'i„,„2\\£{F)\ 



^^^^ " (1 + 2.)i^(«)i V f^illM 

V y(2 + y) i + 2y y 



l + 2y 

2 



□ 



By evaluating equation |H1 at x = y = 1, we derive a new evaluation of the Tutte polynomial 
that counts the number of subtraffics K on G. 

Corollary 5.7. The number of subtraffics on G is equal to S^^^^^hd^, |). 
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